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ABSTRACT 

Current surveys are underway to utilize gravitational lensing by galaxy clusters 
with Einstein radii > 35" in the search for the highest redshift galaxies. Associated 
supernova from the epoch of reionization would have their fluxes boosted above the 
detection threshold, extending their duration of visibility. We predict that the James 
Webb Space Telescope (JWST) will be able to discover lensed core-collapse supernovae 
at redshifts exceeding z = 7-8. 
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1 INTRODUCTION 

Clusters of galaxies act as gravitational lenses, focusing 
light-rays from sources behind them and magnifying their 
images. As this effect enables observers to probe higher red- 
shifts than ever probed before, surveys are being conducted 
with the Hubble Space Telescope (HST) to obtain deep im- 
ages of the sky through massive galaxy clusters. One such 
ongoing program is the Cluster Lensing and Supernova sur- 
vey (CLASH), which is imaging 25 clusters each to a depth 
of 20 orbits (Postman et al. 2012). The 5 clusters selected 
for this program have large Einstein radii of 35" to 55", 
maximizing their potential for discovering ultra-high red- 
shift galaxies. Indeed, three candidate galaxies at redshifts 
z « 9-10 and another candidate galaxy at 2 ~ 11 have al- 
ready been found in the CLASH fields (Bouwens et al. 2012; 
Coe et al. 2013). Similarly, the planned HST Frontier Fields^ 
program will target 6 strong lensing galaxy clusters to reveal 
yet higher redshifts galaxies. 

The James Webb Space Telescope (JWST), the succes- 
sor to HST scheduled for launch in 2018, is likely to have 
analogous observational programs with comparable integra- 
tion times on a similar number of lensing clusters. Although 
the CLASH survey does aim to detect Type la supernova 
(SN) out to redshifts of z ~ 2.5, the current HST cluster 
observations are unlikely to detect gravitationally lensed SN 
from the epoch of reionization at z > 6. Indeed, transient 
science was not identified as a science priority for the Fron- 
tier Fields program, which will not revisit the same field 
twice. The greater sensitivity of JWST and its optimization 
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for observations in the infrared could potentially allow it to 
find lensed supernova from the cosmic dawn in these same 
cluster fields. 

In this Letter, we estimate the cosmic star formation 
rate during the epoch of reionization by requiring that 
enough Pop II stars were formed to ionize the universe. Us- 
ing model spectral time series for Type II SN, as well as a 
simple isothermal sphere model for lensing, we calculate in 
§2-5 the required magnification and duration of detectability 
of such SN at > 6 for different JWST bands and integra- 
tion times. Combining the above, we derive the snapshot 
rate, i.e. the expected number of gravitationally lensed core 
collapse SNe detected in the field-of-view of JWST around 
these high magnification clusters. 



2 STAR FORMATION & SUPERNOVA RATE 

We infer the volumetric supernova rate Rsn{z) as a function 
of redshift by relating it to the cosmic star formation rate 
density (SFRD) p^z): 



RsNiz) = pt{z)rjsN ~ p*iz) 



M„ 



0.7/(, J M V(M) dM 



(1) 



where we use a Salpeter initial mass function (IMF), 
ip{M) oc M~^'^^, and include a factor of 0.7 in the mass 
integral to account for the shallower slope at M < 0.5Mq in 
a realistic IMF (Fukugita, Hogan & Peebles 1998). For the 
stellar mass range between Mmin = 8Mq and Mmax = 40Mq 
appropriate for optically-luminous core-collapse supernova, 
the conversion coefficient between the star formation rate 
and the supernova rate is risN ~ 0.0097Mq^. 
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We require that enough massive stars were formed by 
the end of reionization so as to produce sufficient ionizing 
UV radiation to ionize the intergalactic medium by z^nd = 6. 
This follows the approach used in Pan, Kasen & Loeb (2012), 
albeit with different parameters to bring our estimates closer 
to other inferences in literature, as detailed below. The star 
formation rate during reionization peaks at late times, when 
metals expelled from a prior generation of star formation 
enriched the interstellar gas, so we assume that early Pop 
II stars {Z = O.O2Z0) with a present-day IMF dominated 
the ionizing photon budget. Using the stellar ionizing fluxes 
of Schaerer (2002), we find the average number of ionizing 
photons produced per baryon incorporated into a Pop II 
star was 77^ — 5761. Thus, the mass in stars per comoving 
volume p*{z) should satisfy 

P*{Zend) fj^ fesc = C pb, (2) 

where C is the number of ionizing photons necessary to ion- 
ize each baryon after accounting for recombinations, pb is 
the cosmic baryon density, and fesc is the average escape 
fraction of ionizing photons from their host galaxies into the 
intergalactic medium. Also, we can relate the mass in stars 
per volume p-k{z) to the mass in virialized halos per volume 
via a star formation efficiency 

where we use the Sheth-Tormen mass function of halos for 
dn/dM (Sheth & Tormen 1999), and Afmin ~ IO^Mq is the 
minimum halo mass with atomic hydrogen cooling. The cos- 
mological parameters, such as the matter and baryon den- 
sities Qai, Qb, were taken from Planck Collaboration et al. 
(2013). Assuming is constant, we can calibrate via 
equations (2), (3), and then evaluate p+(2) at any redshift. 
The star formation rate is simply, dpi,{z)/dt. 

Figure 1 shows our estimated SFRD, with C = 3 and 
fesc = 0.2, resulting in a SFRD f» 2 x lO^^M© yr"^ Mpc"^ 
(comoving) between redshifts of 2 = 6 to 8. This corresponds 
to volumetric rates of approximately 2 x 10^* yr^^ Mpc"'^ 
for core-collapse supernova. Our simple SFRD model and 
the resulting SN rates linearly scale with C and /^ci so the 
JWST snapshot rates calculated later can be easily scaled 
for different parameter choices of the SFRD. 



3 LIGHT CURVES 

We adopt the spectral time series of a Type II plateau SN 
from a red giant progenitor with an initial mass I5M0 , com- 
puted by Kasen & Woosley (2009) using a code that solves 
the full multi-wavelength time-dependent radiative transfer 
problem. We plot the SN light curves in the observer frame 
for the best possible HST and JWST filters in Figure 2. 
Note that Type II SN are diverse transients with peak lumi- 
nosities that can vary by more than an order of magnitude, 
and the relationship between the progenitor mass and the 
brightness of the supernova is uncertain; we adopt a single 
characteristic model to represent all core collapse SNe for 
the sake of simplicity. Type IIP SNe are the most common 
events, and the model light curves and spectra used here 
agree very well with observed SNe of average luminosities. 
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Figure 1. Star formation rate density (SFRD) at high redshift. 
The black line shows our fiducial SFRD model used in later cal- 
culations. For comparison, the blue and green regions are taken 
from Robertson & Ellis (2012). The blue region (top) spans the 
high and low values for parametrized star formation histories con- 
sistent with GRB-derived star formation rates, whereas the green 
region (bottom) denotes the SFRD histories derived from UV 
galaxy luminosity densities observed at high redshift, integrated 
down to the observation magnitude limit of Mab ~ —18. Note 
that the latter SFRD is likely to be significantly lower than the 
true cosmic SFRD, as the steep faint-end slope of lower luminosity 
galaxies (possibly down to M^b Sd —10) are omitted (Robertson 
et al. 2013; Ellis et al. 2013), while the GRB-derived SFRD is 
much less flux limited and likely more accurate. Our SFRD pa- 
rameters (C = 3 and fesc = 0.2) were chosen conservatively to 
be consistent with the low end of the GRB-derived SFRD. 

We verified that HST is incapable in practice of de- 
tecting a core-collapse SNe from the epoch of reionization. 
The sensitivity of the HST 1.6^m filter is only a factor of 
2 worse than the JWST F444W filter, but its overwhelming 
drawback is its waveband, which can only probe the SN rest- 
frame UV fiux at 2 ^ 4. Although the JWST F356W is more 
sensitive, the F444W band will be optimal for detecting the 
highest redshift SN that gravitational lensing could provide. 
Figure 3 shows the magnification necessary to detect Type 
II supernova at high redshifts for different integration times. 
Even with a lO^s exposure, a large magnification factor of 
p ^ 10 will be necessary for detecting Type IIP SNe at 
2 > 10 with JWST. 



4 LENSING MAGNIFICATION 

For simplicity, we adopt a singular isothermal sphere (SIS) 
model for the mass distribution of the lensing cluster, within 
which the magnification properties are uniquely specified by 
the Einstein radius 9e (Schneider, Ehlers & Falco 1992). We 
denote the angular separations of the source and the image 
from the center axis of the lens as /3 and 0, respectively. If 
the source lies within the Einstein radius j3 < 6e, two images 
are created at locations 9± = P ± 9e, with magnifications 
p± = 1 ± 9e/ 13. Note that /i_ has negative magnification. 
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Figure 2. Observer frame light curves for a Type IIP supernova from a 15Mq red giant progenitor, for the HST Wide Field Camera-3 
1.6/im filter, and the JWST Near Infrared Camera (NIRCam) F356W and F444W wideband filters at 3.56 and 4.44/jm, respectively. The 
dashed, full, and dotted horizontal lines denote the AB magnitude limits for a IOct detection with 10"*, 10^, and lO^s integration times, 
respectively, for each filter (corresponding to flux limits 50, 13.8, 24.5 njy, respectively, for 10*s exposures). Even a Hubble Deep Field 
measurement has no hope of seeing a regular Type II SN at z ^ 6. A lO^s exposure with JWST can detect a z = 6 supernova without 
magnification. Gravitational lensing would extend its reach to higher redshifts and, more importantly, extend the duration for which the 
supernova remains above the telescope detection threshold. 
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Figure 3. Required magnifications /ir for detecting Type IIP su- 
pcrnovac with JWST at high redshifts. The blue and red lines 
denote the results for the F356W and F444W JWST bands, re- 
spectively, while the dashed and solid lines correspond to integra- 
tion times of 10*8 and lO^s. The latter integration time is similar 
to that used in CLASH. 

that is, the image is flipped compared to the source. If the 
source lies outside the Einstein radius (3 > 6e, there is only 
one image a.t 9 = 9+ with magnification 1 < fx+ < 2. We 
conservatively consider only the higher-magnification image 
at 9+, for which the source angle /3 = 9E/{fJ. — 1). 

Then, the differential source volume (comoving) of mag- 
nified events as a function of magnification and redshift is: 

dV(z, fi) = dA(z, ^i) dDc (4) 

where the differential comoving distance dDc{z) is 
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Figure 4. Comoving source volume as a function of magnification 
II and redshift z over a redshift interval of Az = 1 for a SIS 
lens. The black and green lines denote Einstein radii of 35" and 
55", respectively, while the solid and dashed lines denote z = 6 
and z = 10, respectively. The results are in general agreement 
with more realistic estimates of the search areas per magnification 
factor for the magnification maps of the lensing clusters in the 
CLASH survey (Bouwens et al. 2012). 



with E{z) ~ + z)^ + f^A, and the differential source 
area is 

dA(z,ii) = (2-KDA{z)PDA{z)dp)(l + zf 

= (^27:j-^di?jDA{zf{l + zf. (6) 

Here Da{z) is the angular diameter distance, and the extra 
(1 -I- z)^ is to adjust the area to comoving units. In Figure 4, 
we plot the source volume for a range of Einstein radii typ- 
ical of high-magnification clusters. Given core-collapse SN 
rates of ~ 10"'' yr~^ Mpc""', capturing SN with high mag- 
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nifications within source volumes < 10^ Mpc'' is unlikely. 
Hence, we expect most lensed supernova detected to have 
their fluxes moderately boosted with p < 5; the benefit of 
lensing is to probe somewhat deeper redshifts, and to greatly 
extend the duration of visibility. Also, since high-redshift 
observations are background-limited, for a target signal-to- 
noise ratio, the limiting flux is proportional to t~^^'^, so even 
a modest magnification of p 3 can reduce the required in- 
tegration time by an order-of-magnitude. 

This volume limitation of lensing also justifies our fo- 
cus on core-collapse SNe, which have the highest volumetric 
rates. Although Type la SNe are brighter, their volumetric 
rate is a factor of 4 smaller than the core collapse rate at 
z ~ 7 (Pan, Kasen & Loeb 2012), with the difference dras- 
tically increasing with redshift due to the long delay times 
needed between star formation and explosion for some Type 
la events (Maoz, Mannucci & Brandt 2012). Pair- instability 
SNe from Pop III stars have volumetric rates at least two 
orders of magnitude lower. 



5 SNAPSHOT RATE 

The snapshot 'rate' is the total number of events observed at 
a limiting flux within a given fleld (not per unit time). The 
differential snapshot rate can be calculated from equations 
(1) and (4) via 

N{z, fi) dz dfi = Rsn{z) t{F„, fi, z) dV{z, /i), (7) 

where t{F^,^,z) is the rest-frame duration over which an 
event with magnification jj, will be brighter than the limiting 
fiux Fi, at redshift z, for the observation wavelength v under 
consideration. We find t{Fv,^i,z) using our spectral time- 
series for the Type IIP SN model described in §3. As we 
care about the apparent SN rate for observers, there is an 
implicit factor of (l + z)~^ in front of the intrinsic volumetric 
supernova rate Rsn{z), but that cancels with a (I-I-2) factor 
for t(Fv, jJ,, z) due to cosmic time dilation. 

In Figure 5, we plot the expected snapshot rate of mag- 
nified core-collapse SN detected by JWST above target red- 
shifts, calculated by integrating equation (7) over p and par- 
tially over z. Since NIRCam has two modules each with a 
2.2x2.2 arcmin^ field-of-view, we limit the source area in 
equation (4) to images that lie within this field-of-view. 

We find that a lO'^s JWST snapshot with the F444W 
filter is expected to detect ~ 1 magnified core collapse SN 
at z > 7 around each cluster, and ~ 0.1 SNe at z > 8. 
Using ~ 5 clusters with 0e ^ 35", the prospects for de- 
tecting a few non-superluminous SNe at high redshifts via 
lensing are high. If the other ~ 20 galaxy clusters in the 
CLASH survey with smaller Einstein radii of 8e ~ 15"-30" 
are also included, the expected number of gravitationally 
lensed high-z SNe detected should double. 



6 DISCUSSION 

At z > 6, the observed duration of gravitationally-lensed 
core-collapse SNe can reach > 1 year, lending their detection 
to a search strategy of taking images separated by ~ 0.5 — 1 
year, and looking for flux differences between consecutive 
snapshots. Ideally, the cluster survey should cover most of 




Figure 5. The snapshot rate of gravitationally lensed core col- 
lapse SNe witli JWST, for a single SIS lens with an Einstein 
radius 9^ = 35". Despite the higher sensitivity of the F356W 
band, the F444W band is better for finding lensed SNe at 2 > 6, 
as the SNe remain above the flux limit for a longer time. Note 
that 5-6 high-magnification galaxy clusters with 35" ^ 9e ^ 55" 
are targeted in strong lensing surveys such as CLASH and HST 
Frontier Fields. 

the critical curve area, and not just known locations of mag- 
nified images of high-z galaxies, as the lensed SN may appear 
in currently 'dark' critical curve areas, and serve as a flag for 
its fainter host galaxy. The spectral energy distribution of 
Type II SNe is sufficiently different from blackbody to allow 
for photometric redshift determination, however, typing the 
SNe accurately may require time-consuming spectroscopy. 

Our quantitative results improve upon previous calcu- 
lations of the frequency of lensed SNe. For example, Marri, 
Ferrara & Pozzetti (2000) flrst explored the effects of gravi- 
tational lensing on high-z Type II SNe by intervening cosmo- 
logical mass for different cosmologies, but the predicted de- 
tection rates were unrealistically high because of optimistic 
assumptions about JWST capabilities. Gunnarsson & Goo- 
bar (2003) explored the lensing by massive clusters of distant 
Type la and Type II SNe observed at wavelengths of 0.8- 
1.25 microns, but found the discovery rate tapered off at 
z ~ 3. Also, gravitational lensing is not required per se to 
detect Type II SNe from the epoch of reionization. A mod- 
erate JWST blank-field survey can obtain similar snapshot 
rates, albeit trading off the highest redshift events for more 
lower redshift ones compared to a lensing survey. For exam- 
ple, Mesinger, Johnson & Haiman (2006) found that a lO^s 
exposure with JWST can detect 4-24 SNe per field at z > 5, 
ahhough the assumed SFRD ~ O.IMq yr"^ Mpc"^ was an 
order-of-magnitude higher than our estimates here, and the 
current specifications for JWST NIRCam filter sensitivities 
are now ~ 3 times worse than the values assumed at that 
time. 

For more luminous SNe, Whalen et al. (2012) found 
that core collapse SNe from Pop III progenitors in the earli- 
est galaxies could be visible with the deepest JWST surveys 
(reaching Mab = 32) even at z > 10, as these SNe are 
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bluer and almost an order-of-magnitude brighter than the 
average Type II SNe considered in this paper. Whalen et al. 
(2013) also found that superluminous Type Iln SNe pow- 
ered by circumstellar interactions from Pop III stars could 
be visible out to 2; ~ 20. Truly massive Pop III stars with 
masses M > 200Mq cam also die as extremely bright pair- 
instability supernova, detectable with JWST at 2; > 15 (Pan, 
Kasen & Loeb 2012; Hummel et al. 2012); indeed, the cur- 
rent record for the highest-redshift supernova over observed 
is likely a pair-instability or pulsational pair-instability event 
aX z = 3.90 (Cooke et al. 2012). However, the small volumet- 
ric density of Pop III stars makes it unlikely that these events 
will be strongly lensed. Finally, there is growing evidence of a 
prompt population of Type la SNe, so their volumetric rates 
during the later stages of reionization may not be negligible. 
With the fiducial SFRD model in this Letter, we estimate 
> 1 gravitationally lensed Type la SNe could be discovered 
at 2; > 7 in the snapshots across the ~ 5 high-magnification 
clusters at any given time. 

At lower redshifts, the measured core collapse SN rate 
is a factor of ~ 2 lower than that predicted from the cosmic 
star formation rate (Horiuchi et al. 2011); the most likely 
explanation is that some SN are dim, whether intrinsically 
faint or due to dust obscuration. This will reduce our pre- 
dicted snapshot rate. However, we ignored the contribution 
of multiple lensing images in our analysis. Due to the grav- 
itational lens time delay, which could be ~ 1 — 100 years 
for strong lensing around the clusters of interest (Coe et al. 
2013), multiple images arriving at different times can in- 
crease the expected snapshot detection rate of separate SN 
within the same field-of-view. Although our SIS lens model 
can produce a maximum of only 2 magnified images, sub- 
structure and ellipticity in actual galaxy clusters will likely 
increase both the number of images and their magnifications. 
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